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Abstract 

A dedicated test of the perturbative QCD NLO parton evolution in the very small-x 
region is performed. We find a good agreement with recent precision HERA-data 
for F2(x,Q^), as weU as with the present determination of the curvature of 
Characteristically, perturbative QCD evolutions result in a positive curvature which 
increases as x decreases. Future precision measurements in the very small x-region, 
X < 10~^, could provide a sensitive test of the range of validity of perturbative 
QCD. 



Parton distributions /(x, Q^), f = Q, 9, underlie a Q^-evolution dictated by per- 
turbative QCD at ^ 1 GeV^. It was recently stated 1] that the NLO perturbative 
QCD Q^-evolution disagrees with HERA data [21 El on F'2{x, Q"^) in the small-x region, 
X ^ 10~^. In view of the importance of this statement we perform here an independent 
study of this issue. In contrast to we shall undertake this analysis in the standard frame- 
work where one sets up input distributions at some low Ql, here taken to be Ql = 1.5 
GeV^, corresponding to the lowest considered in [T], and adapting these distributions 
to the data considered. In the present case the data considered will be restricted to 

1.5 GeV^ <Q^ <12 GeV^ 3 x 10"^ <x<3x 10"^ (1) 

as in T] and will be taken from the corresponding measured -^^(a;, Q'^) of the HI 
collaboration Pj. The choice of these data is motivated by their higher precision as 
compared to corresponding data of the ZEUS collaboration 3 , in particular in the very 
small-x region. 

We shall choose two sets of input distributions based on the GRV98 parton distribu- 
tions PJ. In the first set we shall adopt u^, d^, s = s and A = d — u from GRV98 and 
modify u + d and the gluon distribution in the small-x region to obtain an optimal fit to 
the HI data |2] in the aforementioned kinematical region. We shall refer to this fit as the 
'best fit'. The second choice will be constrained to modify the GRV98 m + (i and g distri- 
butions in the small-x region as little as possible. We shall refer to this fit as GRVmod- 
It will turn out that both input distributions are compatible with the data to practically 
the same extent, i.e. yielding comparable x^/dof. In view of these observations we do not 
agree with the conclusions of ref. i.e. we do not confirm a disagreement between the 
NLO Q^-evolution of f{x,Q'^) and the measured [210] Q^-dependence of F^{x,Q'^). 

The remaining flavor-singlet input distributions at Qq = 1.5 GeV^ to be adapted to 
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the recent small-x data are expressed as 

xg{x, Ql) = Ngx-"'' (1 + Agy^ + 7.283x) (1 - x)'-^^^ (2) 
x{u + d) (x, Ql) = Nsa-"^ (l + A^y/x - 4.046a;) (1 - a;)^-^^^ (3) 

where the parameters relevant for the large x-region, x > 10~^, which is of no relevance 
for the present small-x studies, are kept unchanged and are taken from, e.g. GRV98 j3]. 
The refitted relevant small-x parameters turn out to be 

'bestfit' : Ng = 1.70 , 0^ = 0.027, = -1.034 

A^, = 0.171, a, = 0.177, A, = 2.613 (4) 

GRVinod : Ng = 1.443 , ag = 0.125 , Ag = -2.656 

A^, = 0.270 , as = 0.117 , A, = 1.70 (5) 

to be compared with the original GRV98 parameters [4 : A^^ = 1.443, ag = 0.147, 
Ag = —2.656 and A^^ = 0.273, = 0.121, As = 1.80. The resulting predictions are 
compared to the Hl-data [2] in Fig. 1. These results are also consistent with the ZEUS- 
data ini with partly lower statistics. The corresponding x^/dof are 0.50 for the 'best fit' 
{dof = 48) and 0.94 for GRVmod {dof = 50), respectively. Our treatment of the heavy 
flavor contributions to F2 differs from that in We evaluate these contributions in the 
fixed flavor / = 3 scheme of |4^, together with the massive heavy quark (c, b) contri- 
butions, rather than in the f = 4: (massless) scheme utilized in [T- We have checked, 
however, that our disagreement with POj does not result from our / = 3 plus heavy quarks 
vs. the / = 4 massless quark calculations in pQ: we have also performed a fit for / = 4 
massless quarks and the results for F2 and its curvature, to be discussed below, remain 
essentially unchanged. 



In Figs. 2 and 3 we show our gluon and sea input distributions in (2) and (3), as well 
as their evolved shapes at = 4.5 GeV^ in the small-x region. It can be seen that both 
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of our new small-x gluon distributions at = 4.5 GeV^ conform to the rising shape 
obtained in most available analyses published so far, in contrast to the valence-like shape 
obtained in [l^ where the gluon density xg decreases as x — > 0. It is possible to conceive 
a valence-like gluon at some very-low scale, as in but even in this extreme case 
the gluon ends up as non valence-like at Q"^ > 1 GeV^, in particular at = 4.5 GeV^, 
as physically expected. 

Turning now to the curvature test of F2 advocated and discussed in P, we first present 
in Fig. 4 our results for F2{x, Q"^) at x = 10~^, together with two representative expecta- 
tions of global fits [3 E] , as a function of P 

This variable has the advantage that most measurements lie along a straight line as 
indicated by the dotted line at x = 10~^ in Fig. 4. The MRSTOl parametrization j3] 
results in a sizable curvature for F2 in contrast to all other fits shown in Fig. 4. This large 
curvature, incompatible with the data presented in [T], is mainly caused by the valence- 
like input gluon distribution of MRSTOl at Qq = 1 GeV^ in the small-x region which 
becomes even negative for x < 10~^ [Sj. A similar result was obtained in [1^ based on a 
particular gluon distribution xg{x,Q'^) which decreases with decreasing x for x ^ 10^'^ 
even at = 4.5 GeV^ (cf. fig. 7 in pP). More explicitly the curvature can be directly 
extracted from 

F2{x, Q"^) = ao(x) + ai{x)q + a2{x)q^ . (7) 

The curvature a2{x) = ^dgF2{x,Q'^) is evaluated by fitting the predictions for F2{x,Q^) 
at fixed values of x to a (kinematically) given interval of q. In Fig. 5(a) we present 02 (x) 
which results from experimentally selected g-intervals P]: 

0.7 < g < 1.4 for 2 x 10"^ < x < 10^^ 
0.7 < g < 1.2 for 5 x 10"^ < x < 2 x 10"^ 

0.6 < g < 0.8 for x = 5 x 10"^ . (8) 
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Notice that the average value of q decreases with decreasing x due to the kinematically 
more restricted range accessible experimentally. For comparison we also show in 
Fig. 5(b) the curvature a2{x) for an a;-independent fixed g-interval 

0.6 < g < 1.4 (1.5 GeV^ < < 12 GeV^) . (9) 

Apart from the rather large values of a2{x) specific for the MRSTOl fit as discussed 
above (cf. fig. 4), our 'best fit' and GRVmod results, based on the inputs in (4) and (5), 
respectively, do agree well with the experimental curvatures as calculated and presented 
in pP using HI data. It should be noted that perturbative NLO evolutions result in a 
positive curvature a2{x) which increases as x decreases. This feature is supported by the 
data shown in fig. 5(a); since the data point at x < 10~^ is statistically insignificant, 
future precision measurements in this very small x-region should provide a sensitive test 
of the range of validity of perturbative QCD evolutions. 

Furthermore, the HI collaboration [?] has found a good agreement between the per- 
turbative NLO evolution and the slope of -^2(0;, Q^), i.e. the first derivative c?q2F2. 

To conclude, the perturbative NLO evolution of parton distributions in the small-x 
region is compatible with recent high-statistics measurements of the Q ^-dependence of 
F|'(a;, Q^) in that region. A characteristic feature of perturbative QCD evolutions is a 
positive curvature 02 (x) which increases as x decreases (cf. fig. 5). Although present 
data are indicative for such a behavior, they are statistically insignificant for x < 10^^. 
Future precision measurements and the ensuing improvements of the determination of 
the curvature in the very small x-region should provide further information concerning 
the detailed shapes of the gluon and sea distributions, and perhaps may even provide a 
sensitive test of the range of validity of perturbative QCD. 

This work has been supported in part by the 'Bundesministerium fiir Bildung und 
Forschung', Berlin/Bonn. 
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Figure 1: Comparison of our 'best fit' and GRVmod results based on (4) and (5), respec- 
tively, with the HI data 0. To ease the graphical representation, the results and data 
for the lowest bins in = 1.5 GeV^ and 2 GeV^ have been multiplied by 0.75 and 0.85, 
respectively, as indicated 
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Figure 2: The gluon distributions at the input scale Ql = 1.5 GeV^ corresponding to (2) 
with the 'best fit' and GRVmod parameters in (4) and (5), respectively, and at = 4.5 
GeV^. For comparison, the original GRV98 results |1] are shown as well by the dotted 
curves 
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Figure 3: The sea distribution x{u + d) at the input scale Ql = 1.5 GeV^ in (3) with the 
'best fit' and GRVmod parameters in (4) and (5), respectively, and at = 4.5 GeV^. For 
comparison, the original GRV98 results j3] are shown as well by the dotted curves 
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Figure 4: Predictions for F2(x, Q^) at x = 10~^ plotted versus g defined in (6). Repre- 
sentative global fit results are taken from MRSTOl j3] and CTEQ6M jH]. Most small-x 
measurements lie along the straight (dotted) line jT] 
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Figure 5: The curvature 02(3;) as defined in (7) for (a) tlie variable g-intervals in (8) and 
(b) tlie fixed g-interval in (9). Also shown are the corresponding MRSTOl resuhs |5]. 
The experimental curvatures (squares) shown in (a) are taken from 
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